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Abstract 

We consider a problem of shuffling a deck of cards with ordered labels. Namely we split the 
deck of N = fc*q cards (where t > 1 is maximal) into k equally sized stacks and then take the 
top card off of each stack and sort them by the order of their labels and add them to the shuffled 
stack. We show how to find stacks of cards invariant and periodic under the shuffling. We also 
show when gcd(q, k) = 1 the possible periods of this shuffling are all divisors of orders (TV — q). 

1 Introduction 

There are many ways to shuffle a deck of cards. One of the most common is to split the deck into 
equal halves and then riffle or dovetail shuffle the cards back together, wherein the cards from the 
two halves interlace. A perfect shuffle of this type is one where the cards alternate perfectly between 
the two halves. There are two types of perfect shuffles, depending on which card ends up on top. 
These are called "in" and "out" shuffles and have been frequently described in the literature (see, 
e.g. [J or [2]). 

Another way to think about this shuffling is we split the deck in two equal stacks and go through 
the stack from top to bottom using a rule about how to put the two cards into the newly shuffled 
stack. In and out shuffles correspond to the two rules where either we always put the card from the 
second half of the stack on top or we always put the card from the first half of the stack on top. 

In this paper we start the examination of a new kind of shuffling where our rule is to set an 
order on the labeling of the cards (we allow for labels to be used multiple times in the deck) and 
we again go through from top to bottom but now we let the order of the labeling on the cards 
determine which one goes on top. That is, starting with N = kn cards with labels from the ordered 
list A\ y A2 >- • • • >- Aj, divide the stack of cards into k stacks of n, then take the top card off of 
each stack, sort the k cards (according to the order of the labels, where if the labels agree then the 
order they are sorted is unimportant) and add them to the new stack. An example of this is shown 
in Figure [T] where the labels are 2 >- 1 y and we let the heaviest weights "sink down" (i.e., go to 
the lower card height where we start counting from the top card down). 

Equivalently, this is the same as starting out with a list of kn labeled objects, 

ao, ai, . . . , a n _i, a n , a n+ i, . . . , am-\, a2n, ■ ■ ■ , «fcn-i, 
and putting this into a kxn matrix where we proceed by filling up the rows left to right and top 
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Figure 1: An example of shuffling with ordered cards taking the deck 01211201 to the deck 10212011. 



to bottom, 



Ol 



a. 



«2n-l 



\ «fen-l/ 

Now take this matrix and sort the elements in each column according to the ordering of the labels, 

.. . \ 



/ b 



bk 
bk+i 



\bk-i fok-i ■ ■ ■ bkn-ij 
and finally concatenate the columns to form the new list. 

bo,h,. ■ ■ , bk-\,bk, bk+i, • • • , b 2 k-i,b2k, ■ ■ ■ , bkn-i- 

As an example if we have k = 3 and start with iV = 12 cards labeled 021100122110 (where we 
again let the 3 labels be ordered 2 >- 1 >~ 0) then we have 



021100122110 

This process can be repeated and so we have 
021100122110 — ► 200210111210 




200210111210. 



211200110210 



200210111210 



and now we see we have a stack which will repeat itself every two shuffling steps, or in other words 
a periodic stack. A stack which returns to itself after one shuffling step will similarly be called a 
fixed stack. 

In general, we can relate the shuffling to a directed graph where each of the possible j N stacks 
are the vertices and we put a directed edge between two vertices if shuffling one stack gives the 
other. Since the outcome of our shuffling is uniquely determined by the order of the cards, each edge 
will only have one edge going out (though it is possible many edges can go in). This immediately 
gives us the following. 

Observation 1. Given any stack of cards after applying the shuffling procedure finitely many times 
we will settle into a stack which is periodic under the shuffling. 
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There now arise several natural questions. For example, how do we find fixed/periodic stacks? 
What periods are possible? How many periodic/fixed stacks are there? How long does it take for 
a stack to settle into a periodic stack? 

In this paper we will answer some of these questions. In Section [2] we will introduce a weight 
function on the subscripts and show how to use this to represent the shuffling by a poset structure. 
In Section [3] we will show in the case when gcd(q, k) = 1 (where N = k l q) the possible periods are 
all divisors of order&(iV — q). In Section [4] we will show how to adopt the poset structures to find 
posets that generate all fixed and periodic stacks. Finally, in Section [5] we will give some concluding 
remarks. 

We will throughout assume the number of cards is N = k l q = kn, where t > 1 is the highest 
power of k that divides N, and n = k t ~ 1 q is the size of the stacks we split N into when shuffling. 
For simplicity we will focus on subscripts, i.e., i —> j means — > bj. 



2 Representing our shuffling using a poset structure 

The key to understanding this shuffling is looking at a column which will (for some £ € {0, 
consist of the terms 

£,£ + n,£ + 2n,...,£ + (k-l)n, 

and after sorting will be sent to 

k£,k£ + l,k£ + 2, . . . ,k£ + (k - 1). 
For example, returning to the case N = 12 and k = 3 then we get the following. 



,71-1}) 




4~5] |6 7 8| 9 10 111 



This shows, for example, {00,04, as} — > {^o^i;^} in some order (depending on the labels of the 
cards). To help us understand what is happening it is useful to weight the subscripts. 

Definition 1. A shuffling weight function on the subscripts is a map 92 : {0, . . . , N — 1} — > 7L which 
satisfies the following two conditions for £ 6 {0, 1, . . . , n — 1}; 

(i) {<p(£), <p(l + n), . . . , ip(£ + (k - l)n)} = {(p(k£), ip{kl + 1), . . . , ip(k£ + (k - 1))}. 

(ii) <p(k£) < ip(M + 1) < • • • < <p(k£ + (k— 1)). 

The first condition says the weight of the entries in the column and the weight of the entries in 
a block it maps to are equal. The second condition says the weights are distinct and increasing in 
a block (and combined with the first condition says the weights of the column are also distinct). 

As an example for N = 12 and k = 3 one possible weight function is given below. 
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If we now take this weight function and replace all the entries in the above diagram with the 
corresponding weights then our diagram now becomes the following (for which it is easy to check 
the two conditions are satisfied). 




12| 


10 1 2| 


10 1 2| 


1 2 



We can use this weight function to map the subscripts to the subscripts in a bijective manner 
so the weight is preserved. For instance in the third column we have {2,6,10} — > {6,7,8}. Since 
^(2) = 2 = ip(8), ip(6) = = ip(6) and ^(10) = 1 = ip(7) then we would have 2 -»• 8, 6 -»• 6 and 
10 —¥ 7. Repeating this for each column/block combination we can now break up the elements of 
{0, 1, . . . , N — 1} into cycles, and we can place these cycles into a poset structure where the height 
of a cycle is the weight of a subscript in the cycle (by construction they are all equal). So for our 
particular weight function we would end up with three levels with the following cycles. 



P ^7—10^ 
ffP 4 -9-^3 JI 

This diagram essentially represents half of the shuffling (i.e., how we move from columns to 
blocks in the concatenation). We also need to represent the other half of the shuffling, which is 
sorting among the columns. To do this we add edges between levels in the diagram connecting any 
two subscripts which appear in the same column. Doing this gives us the following diagram (which 
we will refer to as the shuffling poset). 




The important part about the shuffling poset is that the higher ordered cards (or by analogy 
heaviest cards) will always try to "sink down", i.e., the highest label in a column will want to 
map to the smallest possible ip value which corresponds to the lowest level it can reach, the second 
highest label will similarly go to the second lowest level it can reach and so on. In particular, we 
have the following. 

Observation 2. Shuffling can be carried out completely in the shuffling poset by placing the cor- 
responding cards in their position and carrying out the following two steps: 

(i) Using only the edges between levels in the poset swap cards so no edge connects a higher 
ordered card at height (or weight) i in the poset with a lower ordered card at height (or 
weight) j where i > j. 

(ii) Using only the edges in the levels advance the card to the next entry. 
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2.1 Generating a shuffling weight function 

We now see that once we have a shuffling weight function we can generate the shuffling poset which 
completely describes the shuffling. We now show a shuffling weight function always exists by giving 
an algorithm which will generate one. 

(1) Set weight = and construct a kxn array where the element in the (i, j)th entry of the 
array (0 < i < k — 1 and < j < n — l) is a : b = (i + kj) : ((i + kj)%n), where s%t is the 
remainder when s is divided by t. 

(2) While there are still cells in the array which have not been crossed off do the following: 

(i) Construct a directed graph on the vertices {0, 1, ... ,n — 1} by letting j — > k where j is 
the column and k is the b entry in the lowest cell of column j which has not been crossed 
out. (If all the cells in column j have been crossed out then j will be an isolated node.) 

(ii) Find all directed cycles in the graph that was constructed. For each edge in the directed 
cycle let = weight, where a comes from the cell that generated the edge. 

(iii) Cross out all cells which were used to assign a weight. 

(iv) Increase weight by 1. 

An example of the algorithm being carried out when N = 32 and k = 4 is shown in the 
Appendix. (The reader is encouraged to try this algorithm out to generate the weight function for 
N = 12 and k = 3 given above and for N = 24 and k = 6 given in Section [5j) 

Theorem 1. The above algorithm generates a shuffling weight function. 

Proof. First let us observe that the algorithm will terminate so (p is a well defined function. This 
follows since each column is referred to exactly k times and so if all of the entries in a column have 
been crossed out no edge will be directed towards it. In particular, if we ignore isolated vertices then 
the directed graph generated in step 2i will have all vertices with outdegree 1 and so it must contain 
at least one directed cycle (in fact one directed cycle for each connected component). Therefore at 
every stage we will continue to cross out cells unless all cells have been crossed out so the algorithm 
will terminate. 

The first condition of a shuffling weight function (that weights in the blocks match weights in 
the columns) is satisfied because we are pulling out cycles. Namely, if we look at a cell a : b in the 
jth column, this can be understood as the subscript a will be placed into the bth column (when 
shuffling) . What we need is to make sure there is some subscript which will go into the j th column 
that also will be assigned the same weight. But this happens because we assign the entire directed 
cycle the same weight and by step 2ii the vertex preceding j will give a subscript that maps into 
the bth column. 

The second condition of a shuffling weight function (that weights in each block are increasing) 
is easily satisfied since the blocks form the columns of the array generated in step 1, and in columns 
the weights assigned to the cells are increasing (since in each round at most one cell in each column 
will be assigned the current weight.) □ 

The algorithm used the lowest cell which had not been crossed out. We could also have used 
the highest cell which was not crossed out and then decreased the weight by 1. This will essentially 
generate the same picture (where we rotated the arrays by 180 and let i —¥ N — 1 — i in the entries 
and in the directed graphs generated in step 2i of the algorithm. 
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Another important thing to note is the directed cycles we are pulling out are the same directed 
cycles (but with the b terms replaced by the a terms from the corresponding cells) found in the 
shuffling poset. If we combine this with the previous idea of going from bottom to top we get the 
following observation. 

Observation 3. In the shuffling poset i — >■ j if and only ifN — \ — i—>N — 1— j. 

Something else worth noting is that if we look at the shuffling weight function produced in the 
Appendix we see there is some symmetry, i.e., <p{i) + <p(31 — i) is independent of i. So for example 
we have ip(2) + </?(29) = 7 = (p(l2) + <^(19). When this happens we will call this a symmetrical 
weight function. 

Conjecture 1. The shuffling weight function produced by the algorithm is symmetric. 

It is easy to show a symmetric weight function always exists for every N and k. For instance 
let Pup be the weight function generated by the given algorithm (where we use the lowest cells and 
work our way up) and let <~Pdown 

be the weight function generated by the modified algorithm (where 
we use the highest cells and work our way down). Then the weight function ip = (p up + (fdown is 
still a shuffling weight function and by the symmetry of the two algorithms will be symmetric. 

3 A simple weight function when gcd(g, k) = 1 

In the preceding section we saw that a weight function could be used to give us the shuffling poset, 
which in turn gives us another representation of how to shuffle. One natural question to ask is 
what can we say about the lengths of the directed cycles in the shuffling poset. In this section we 
will show in the special case when N = k l q and gcd(c/, k) = 1 this question has an easy answer. 
Note this will cover all the cases when k is a prime (in particular k = 2). First, for gcd(q,k) = 1 
we observe one shuffling weight function can be found using the base k expansion of subscripts. 

Lemma 2. Let N = kn = k l q with gcd(q, k) = 1 and let . . . A t A t -\ . . . A\Aq be the base k expansion 
of A. Then ip(A) = Aq + ■ ■ ■ + A t -\ is a shuffling weight function. 

Proof. Let £ € {0, 1, . . . , n — 1} and . . . LjL^i . . . L\Lq be the base k expansion of £, and let the base 
k expansion of n be ... uO ... 00 where u ^ is in the (t — l)th slot and gcd(-u, k) = gcd(g, k) = 1 
(this follows from n = ^^q). 

In particular, for i G {0, . . . , k — 1} the base k expansion of k£ + i will be 

• • • Lt-iLt-2 ■ ■ ■ Lqi 

which has ip(k£ + i) = i + Lq + • • • + L t -2- On the other hand the base k expansion of £ + in will be 

...((L t -i + iu)%k)L t -. 2 ...L 1 L 

which has ip(£ + iu) = ((L t _i + iu)Vok) + Lq + • • • + L t _2 (where s%t is the remainder of s when 
divided by t). Because gcd(u, k) = I then (L t _i + iu)%k will cycle through all k possibilities as % 
goes from to k — 1. It follows the first condition of a shuffling weight function is satisfied. 

The second condition of a shuffling weight function is satisfied since the base k expansion of 
k£ + i for i G {0, . . . , k — 1} is . . . L t -\L t -2 ■ ■ ■ Loi and so 

ip(kl + i) = i + L + ---+L t - 2 - 
From this it follows ip(k£) < ip(k£ + !)<■■■< (p(k£ + (k — 1)). □ 
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One nice feature about this case is the the rule for mapping subscripts is easy to describe. 
Namely, if the base k expansion of A is . . . At-iAt-2 . . . A\Aq then the map is 

A->kA + At-! (modiV). (1) 

(There are two things to check, one is that both these terms have the same weight and the other 
is that the column that contains A will map to the block that contains kA + At-±, both conditions 
are easily checked.) 

We will now use this map to determine the cycle lengths in the shuffling poset. 

Theorem 3. Let N = k l q with gcd(k,q) = 1, and let orderfc(s) denote the multiplicative order of 
k modulo s. Then the length of a cycle in the shuffling poset when we divide N into k equal stacks 
for shuffling is a divisor of order/% (TV — q) . Further, there is a cycle of length order/% (N — q) in the 
shuffling poset. 

Before we begin we note by our assumption that gcd(iV — k,q) = 1 and so the multiplicative 
order is well defined. As a check we note order3(12 — 4) = order3(8) = 2 since 3 2 = 9 = 1 (mod 8), 
this agrees with the diagram given above for TV = 12 and k = 3. 

Proof. Consider a cycle starting at x, and suppose that the base k expansion of x is . . . A t -\ . . . A\Aq. 
Using ([!]) we have after t steps we will be at 

x ->■ kx + A t -i (mod N) 

-> k 2 x + kA t -x + M-i (mod N) 

-> k 3 x + k 2 A t _ x + kA t -2 + A t -z (mod N) 

-¥ ■■■ 

t-l 

-> fcW^A^A (modiV). 

i=0 




In particular, after we have taken t steps, the last t terms in the base k expansion will agree with 
the last t terms in the base k expansion of x. Now suppose we repeat this r times (so a total of rt 
steps). Then we have 

x -> k l x + A' (mod N) 

k 2t x + k f A' + A' (mod N) 
-> k 3t x + k 2t A' + k*A' + A' (mod N) 
-¥ ■■■ 

r-1 

-> k rt x + ^k u A' (modiV). 

i=0 

For some r we will be back where we started if 

r-1 

k rt x + k u A' = x (mod N = k l q). 
i=0 

Multiplying both sides by k l — 1 and simplifying this is equivalent to 

(k ri - 1) (x(k* - 1) + A') = (mod (fc 1 - l)k*q) 
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Looking at the base k expansion of x we have x = A' + mk 1 for some m, if we now substitute this 
in and simplify we get 

(k rt - l)k* (A' + m(jfe* - 1)) = (mod (Jfc* - l)jfc*qr) 

or 

(k rt - 1) (A' + m(fc* - 1)) = (mod (jfe* - l)g). (2) 

Since (fc* — l)g = N — q, ii rt = orders (iV — g) then fc r< — 1 = (mod (fc* — l)g) and this condition 
is satisfied. In particular, after taking orderfc(iV — q) steps then x — > x for each value of x, and so 
each cycle must be some divisor of orders (N — q). 

We are implicitly using 1 1 orders [N — q). To see why this is true, let M = orders (N — q) then 
N -q=(k t -l)q\ \k M - 1), and so 

(jfe* - 1) | (k M - 1) so gcd(fc* - 1, k M - 1) = k l - 1, 

but we also have in general 

gcd(A; a - 1, k b - 1) = fcS cd (°> 6 ) - 1 so gpd{k* - 1, k m - 1) = yfc gcd(t ' m) - 1. 

Combining these two statements we have gcd(t, M) = t showing t is a divisor of M = orders (N — q) . 

Finally, to show there is a cycle of length order^(A^ — q) we note ([2]) also holds for x = 1, which 
corresponds to A' = 1 and m = 0. So ([2]) reduces to finding the smallest value of rt so 

k rt = l (mod JV-g), 

which is clearly rt = orders (iV — g) . □ 
3.1 A more generalized weight function 

The preceding weight function relied on having gcd(g, k) = 1 (this was used in the lemma to make 
sure we hit all of the residue classes modulo k, and then in the Theorem [3] by giving a simple rule 
for mapping). We now present a weight function that works in more cases (but from the definition 
it will agree with the previous weight function when gcd(g, k) = 1). For example, the following 
weight function will work for the case when k is any prime power. 

Lemma 4. Let N = kn = k l q with gcd (q/ gcd(g, k), gcd(g, k)) = 1 and let . . . A t A t -iA t -2 ■ ■ ■ AiAq 

be the base k expansion of A. Then (p(A) = Ao H h At-\ + [At% gcd(/c, g)) is a shuffling weight 

function. 

Proof. For I € {0, . . . , n — 1} let ... LfLt-i . . . L±Lq be the base k expansion of t. Then for i € 
{0, . . . , k — 1}, the base k expansion of kl + i is ... Lt~\Lt~2 ■ ■ ■ L\L^i and so 

<f(k£ + i) = i + L Q + --- + L t _ 2 + {Lt-i% gcd(k, q)) . 

Using this, we note the second condition for a shuffling weight function is easily satisfied. 

Since n = k l ~ l q, the base k expansion of n is ... NfNt-iO ... 0. Also we have Nt-i = (q%k), 
and so gcd(g, k) = gcd ((q%k), k) = gcd(N t -i, k). So for i £ {0, . . . , k — 1} the base k expansion of 
£ + in is ... CfCt-iLt-2 ■ ■ ■ Lq, i.e., it agrees in the first t — 2 slots with the expansion of £, so we 
need to understand Ct~i + (Ct% gcd(g, k)) . 

For a given value of i we have Ct-\ = ((Lt-i + iN t -i)%k). Since Nt-i/ gcd(q,k) is relatively 
prime to k, this takes on k/ gcd(q,k) different values that differ by a multiple of gcd(g, k) from 
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At-i (modulo k), and we will attain each one of these values gcd(q,k) times as u ranges over its 
possible values. In particular we have, 



C t -i € |(L t _i%gcd(g,fe)),(L t _i%gcd(g,fc)) + gcd(q,k), 

. . . , (L t ^%gcd(q,k)) + (^^y " l) gcd(g,fc)} 

Fix a possible value for Ct—\ and let i be the first value of i that gives us Ct-i, and let Ct be the 
value of Ct for this corresponding i. In particular, the only values of i where we will be at the fixed 
value of Ct-i are i = i + rakj gcd(g, k) for m £ {0, ... , gcd(g, k) — 1}. For these values of u we have 

(The last step follows from noting q = Nt—i + kNt + k Nt+l H and that when we divide this by 

gcd(g, k) all but the first two terms will have at least one factor of k.) So we have 



(C t %gcd(q,k)) = ((^ + m gcd ( g g k) ^%gcd(q,k) 



Since gcd [q/ gcd(q,k),gcd(q,k)) = 1 then this covers all of the residue classes modulo gcd(q,k). 
Combined with what we know about Ct-i then we have 

C t ^+(C t %gcd(q,k)) = 

{ (Lt-i% gcd(g, k)) , (Lt-i% gcd(g, k)) + 1, . . . , (L t ^% gcd(g, k))+k-l}. 

So we have 

tp(£ + in) = j + L + • • • + L t -2 + gcd(/c, g)) , 

where j which ranges over {0, . . . , k — 1} as i ranges over {0, . . . , k — 1}. □ 

The problem is even though we have a simple weight function, the rule for mapping is not 
simple, and so we have no similar result as Theorem [3j We note one of the results coming out of 
the proof of Theorem [3] is when gcd(g, k) = 1 the cycle in the shuffle poset which contains 1 has 
maximal length. This no longer needs to hold when gcd(q, k) ^ 1. If we look at the example shown 
in the Appendix, we can see the cycle which will contain 1 has length 3 but the maximal length of 
a cycle is 6. 



4 Finding fixed and periodic stacks 

In Section [2] we saw a way to represent our shuffling in terms of a poset of cycles. We will now 
exploit this poset to find stacks of cards which are fixed or periodic under shuffling. 

To find the fixed stacks, we observe when looking at the cycles of the shuffling poset in a stack 
which is fixed under shuffling there are two necessary conditions: 

(i) If there is an edge between levels in the shuffling poset then the label on the lower level must 
be at least as great as the label on the higher level. (Otherwise we would swap labels and we 
would not be fixed.) 
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(ii) All the labels in a cycle on a level in the shuffling poset must agree. (Otherwise when we 
shift by one in the cycle the stack is not fixed.) 

It is easy to see that these conditions are also sufficient. 

We now form a poset, which we call the fixed poset, based off the shuffling poset. Namely, 
each cycle goes to one element in the poset (at the same height as before) and we connect an edge 
between two corresponding cycles if one cycle contains c while the other contains d from the same 
column and there is no e in the column with (p(c) < (p(e) < (p(d). (Technically, by condition (i) 
above we would want to connect all cycles which are connected by an edge in the shuffle poset, 
but by transitivity we only need to consider edges which cannot be broken down further.) Some 
examples of fixed posets are given in Figure [2j 




Figure 2: Fixed posets for various values of N and k. 



Using the two conditions, the fixed stacks are now found by assigning labeled cards to each 
element in the fixed poset so the labels are weakly increasing with respect to the fixed poset. This 
allows us to quickly and easily find fixed stacks for a given N and k. 

It also helps us to answer whether or not fixed stacks with some given property can exist. For 
example if N = 52 and k = 2 it is easy to see from the fixed poset there is no fixed stack with four 
labels each label with 13 cards (i.e., such as in a standard deck of cards with the labels being the 
suits in some order). 

Note the fixed posets in Figure [2] are symmetric in that if we flipped it upside down we would 
have the same poset. This will always happen when our weight function is symmetric (and as 
already noted we can always find such a weight function). 

One question to consider is the number of fixed stacks. This is dependent on both the fixed 
poset and the number of labels. As an example, for a fixed k and a fixed order of the labeling we 
can have dramatically different behavior for various values of N. For example, suppose we have 
two labels 1^0 and k = 2. If N = 2*, then by Theorem [3] all the cycles have length at most t and 
there are t + 1 levels in the fixed poset using the weight function given in Lemma [2j We can then 
assign all of the cycles with height < t/2 in the poset to have label 1 and then take any collection 
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of posets from the layer at height \t/2] to have label 1. There are at least 

1 / t \ cA 
t\\t/2\) * (IniV) 3 / 2 ' 

cycles on the \t/2] level for some constant c > 0, and since we can take any subset of them and 
form a fixed stack then there are at least 2 cN /^ n N ) a/2 fixed stacks for A = 2 t (in particular this is 
super-polynomial). By comparison if we consider A = 4-3*, then it is easy to show it has a fixed 
poset of the following form. 

t + 1 cycles 
1 cycle 
t + 1 cycles 

So in particular there are exactly 4-2* < cN 0A05 = o(A) fixed stacks for A = 4-3*. 
To make the comparison more concrete we have for A = 1024 = 2 10 there are 292, 592, 830 fixed 
stacks, while for A = 972 = 4-3 5 there are 128 fixed stacks. 

4.1 Finding periodic stacks 

When looking for periodic stacks we will again set up a poset, but instead of having cycles as 
elements in the poset we will have individual subscripts be the elements. The basic idea is to 
consider what can keep a card at a given subscript from dropping down as we go through repeated 
shufflings (because we are periodic the only operation that will happen as we shuffle is shifting the 
cycles). In particular, an element A can only drop down if for some B and some s > we have the 
following edges in our shuffling poset: 

A = c ->ci-^ >c s -.i^c s 

i 

B = do— >•• ■ ■ — >-d s _i— >d s 

In this case we need to make sure the card in B will not cause the card in A to sink. To do this 
we draw the poset where the elements are the subscripts (as before we can place A at level <p(A)), 
and we connect an edge between A and B if we have the edges in our shuffling poset as indicated 
above and there is no e in the same column as c s and d s so f(c s ) < 93(e) < <p(d s ). We will call this 
the periodic poset. An example of the situation is shown in Figure [3j 

As before, the periodic stacks are now found by assigning labeled cards to each subscript in the 
periodic poset so the labels are weakly increasing with respect to the periodic poset. This allows 
us to quickly and easily find periodic stacks for given and k. 

The possible periods of the periodic stacks are determined by the size of the cycles in the 
shuffling poset. Namely, the possible periods are the divisors of the least common multiple of the 
cycle lengths in the shuffling poset. In the case when gcd(g, k) = 1 then Theorem [3] shows the 
possible periods of stacks are divisors of orders (A — q), on the other hand it is easy to construct a 
periodic stack for any period dividing orders (A — q). 
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(a) Shuffling poset 



(b) Periodic poset 



Figure 3: Posets for the case N = 24 and k = 2. 

5 Concluding remarks 

We have seen how to find a weight function which can in turn allow us to represent our shuffling 
in a shuffling poset. This poset can then be modified to help us find fixed stacks, periodic stacks, 
and also tell us which periods are possible. 

One of the problems we have not addressed is how quickly a stack will settle into a periodic 
stack. As with the number of fixed stacks this depends highly on N. For example, it is not hard to 
see (i.e., using the periodic poset) that it takes no more than 3tm shuffles to settle into a periodic 
orbit (where m is the least common multiple of the cycle lengths and t + 1 the number of levels in 
the shuffling poset). So for example when N = 2* and k = 2 then we need at most 3( lniVJ steps. 
On the other hand for ./V = 4-3 fc it is easy to construct a stack that takes exactly N/2 steps to get 
into a periodic stack. 

There are still many questions that remain. One of the biggest problems is trying to understand 
how the weight function works for arbitrary N and k. For instance the weight functions given in 
Section [3] do not apply for the case N = 24 and k = 6. In this case the algorithm for finding a 
weight function generates the following: 
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One notable difference between this weight function and the weight function when gcd(g, k) = 1 
or gcd (q/ gcd((/, k), gcd(g, k)) = 1 is that in the blocks the weight function does not consist of 
consecutive numbers, i.e., it has gaps. Determining why these gaps are there and where they will 
appear given N and k will go a long way to understanding the shuffling weight function. 

Another important question in regards to the shuffling posets is understanding the possible 
cycle lengths. We understand what is going on for the case when gcd(g, k) = 1, but all other cases 
remain open. For example, is it true the least common multiple of the cycle lengths is the length 
of the longest cycle? 

We can also consider what happens when instead of only considering a single type of shuffling we 
consider combining the j\ different shuffling rules that come from all the possible rearrangements of 
the ordering of the labels. And of course, perhaps the most important thing missing right now is a 
good magic trick that can be performed using this shuffling rule, which was the original motivation 
of Larry Carter and J.-C. Reyes who first suggested this problem! 
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Appendix 

We implement the algorithm given in Section 12.11 to find the weight function in the case ./V 
and k = 4. The steps are shown below. 
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The generated weight function is given in the following table. 
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